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ABSTRACT 

An extensive multi-satellite campaign on NGC 5548 has revealed this archetypal Seyfert-1 galaxy to be in an exceptional state of 
persistent heavy absorption. Our observations taken in 2013-2014 with XMM-Newton, Swift, NuSTAR, INTEGRAL, Chandra, HST 
and two ground-based observatories have together enabled us to establish that this unexpected phenomenon is caused by an outflow¬ 
ing stream of weakly ionised gas (called the obscurer), extending from the vicinity of the accretion disk to the broad-line region. In 
this work we present the details of our campaign and the data obtained by all the observatories. We determine the spectral energy 
distribution of NGC 5548 from near-infrared to hard X-rays by establishing the contribution of various emission and absorption pro¬ 
cesses taking place along our line of sight towards the central engine. We thus uncover the intrinsic emission and produce a broadband 
continuum model for both obscured (average summer 2013 data) and unobscured (< 2011) epochs of NGC 5548. Our results suggest 
that the intrinsic NIR/optical/UV continuum is a single Comptonised component with its higher energy tail creating the ‘soft X-ray 
excess’. This component is compatible with emission from a warm, optically-thick corona as part of the inner accretion disk. We 
then investigate the effects of the continuum on the ionisation balance and thermal stability of photoionised gas for unobscured and 
obscured epochs. 

Key words. X-rays: galaxies - galaxies: active - galaxies: Seyfert - galaxies: individual: NGC 5548 - techniques: spectroscopic 


1. Introduction 

The supermassive black holes (SMBHs) at the heart of active 
galactic nuclei (AGN) grow through accretion of matter from 
their host galaxies. This accretion is accompanied by outflows 
(powerful relativistic jets of plasma and/or winds of ionised 
gas), which transport matter and energy away from the nucleus, 
thus linking the SMBHs to their host galaxies. Understanding 


the feedback mechanisms between SMBHs and their environ¬ 
ments is important in AGN science, as well as in cosmology. 
They can have significant impacts and implications for the 
evolution of SMBHs and star formation in their host galaxies 
(e.g. Silk & Rees 1998; King 2010), chemical enrichment of 
their surrounding intergalactic medium (e.g. Oppenheimer & 
Dave 2006), and the cooling flows at the core of galaxy clusters 
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(e.g. Ciotti & ©striker 2001). 

The ionised outflows are an essential component of the 
energy balance of AGN and are best studied through high- 
resolution X-ray and UV spectroscopy of their absorption line 
spectra. The Reflection Grating Spectrometer (RGS) of XMM- 
Newton and the Low-Energy & High-Energy Transmission 
Grating Spectrometers (LETGS and HETGS) of Chandra, 
together with the Cosmic Origins Spectrograph (COS) of the 
Hubble Space Telescope (HST) have vastly advanced our 
knowledge of these outflows in recent years (see e.g. the review 
by Costantini 2010). With current instrumentation, the nearby 
and bright Seyfert type AGN are the best laboratories for 
studying these outflows. The majority of the absorbing gas is 
detectable in the X-ray band. The outflows detected in the soft 
X-rays (referred to as ‘warm absorbers’) are found to have 
column densities of 10^*^-10^^ cm“^ and consist of multiple 
phases of photoionised gas with temperatures of 10^^-10^ K, 
travelling at velocities of up to a few 10^ km s“' (see e.g. 
Blustin et al. 2005). Understanding the origin and launching 
mechanism of the ionised outflows in AGN is an active area of 
research. It has been suggested that these outflows are produced 
by irradiation of the dusty gas torus structure, which surrounds 
the SMBH and accretion disk (e.g. Krolik & Kriss 2001), or that 
they originate as radiatively driven winds from the accretion 
disk (e.g. Begelman et al. 1983; Proga et al. 2000), or are in 
the form of magnetohydrodynamic (MHD) winds from the 
accretion disk (e.g. Konigl & Kartje 1994; Bottorff et al. 2000). 

Many aspects and physical properties of these outflows 
are still poorly understood. Eor instance, the location of the 
ionised absorbers needs to be established to distinguish between 
the different outflow mechanisms and to determine their mass 
outflow rates and kinetic luminosities, which are essential 
parameters in assessing their impact on their surroundings and 
their contribution to AGN feedback (e.g. Hopkins & Elvis 
2010). The distance of an X-ray or UV absorbing outflow to the 
ionising source can be determined from estimates of the density 
of the absorbing gas. The density can be measured from either 
density-sensitive UV lines (e.g. Gabel et al. 2005) or from the 
recombination timescale of the ionised absorber (Kaastra et al. 
2012). In the X-ray band, the latter delivers more robust density 
measurements than the use of density-sensitive X-ray lines. In 
the recombination timescale method, as the intrinsic luminosity 
of the AGN varies over time, the ionisation state of the absorber 
changes with a time delay; by measuring this lag, the electron 
density and hence the distance of the absorber to the ionising 
source can be obtained. Such a study was first done as part of 
our 2009 multi-wavelength campaign on the Seyfert-1/QSO 
Mrk 509 (Kaastra et al. 2012). Erom the response of three 
absorber components with the highest ionisation in the soft 
X-rays, individual distances of 5-100 pc were derived, pointing 
to an origin in the narrow-line region (NLR) or the AGN torus 
region. 

The archetypal Seyfert-1 galaxy NGC 5548 is one of the 
most widely studied nearby active galaxies. It was one of the 
12 classified objects in the seminal work of Seyfert (1943), and 
since the sixties it has been the target of various AGN studies. 
In more recent times it was the first object in which narrow 
X-ray absorption lines from warm absorbers were discovered by 
Kaastra et al. (2000) using a high-resolution Chandra LETGS 
spectrum. Later on, a detailed study of its warm absorber was 
carried out by Steenbrugge et al. (2005). Eurthermore, extensive 


optical/UV reverberation mapping studies of NGC 5548 (e.g. 
Peterson et al. 2002; Pancoast et al. 2014) have provided one 
of the most detailed pictures available of the broad-line region 
(BLR) size and structure in this AGN. Prior to our recent 
campaign on NGC 5548, the X-ray properties and variability of 
NGC 5548 were known to be typical of standard Seyfert-1 AGN 
sharing common characteristics (e.g. Bianchi et al. 2009; Ponti 
et al. 2012). 

In 2013-2014 we carried out an ambitious multi-wavelength 
campaign on NGC 5548, which was similar but more extensive 
than our campaign on Mrk 509. This remarkable campaign 
has utilised eight observatories to take simultaneous and 
frequent observations of the AGN. It incorporates instruments 
onboard five X-ray observatories: XMM-Newton (Jansen et al. 
2001), Swift (Gehrels et al. 2004), NuSTAR (Harrison et al. 
2013), INTEGRAL (Winkler et al. 2003), Chandra’s LETGS 
(Brinkman et al. 2000), as well as the HST COS (Green et al. 
2012), and two ground-based optical observatories: the Wise 
Observatory (WO) and the Observatorio Cerro Armazones 
(OCA). These observatories have collected over 2.4 Ms of 
X-ray and 800 ks of optical/UV observation time. As previously 
reported by Kaastra et al. (2014), NGC 5548 was discovered 
to be obscured in X-rays with mainly narrow emission features 
imprinted on a heavily absorbed continuum. This obscuration 
is thought to be caused by a stream of clumpy weakly-ionised 
gas located at distances of about 2-7 light days from the black 
hole and partially covering the X-ray source and the BLR. Erom 
its associated broad UV absorption lines detected in HST COS 
spectra, the obscurer is found to be outflowing with velocities of 
up to 5000 km s ^ The intense Swift monitoring on NGC 5548 
shows the obscuration has been continuously present for a 
few years (at least since Eeb 2012). As the ionising UV/X-ray 
radiation is being shielded by the obscurer, new weakly-ionised 
features of UV and X-ray absorber outflows have been detected. 
Compared to normal warm absorber outflows commonly seen 
in Seyfert-Is at pc scale distances, the remarkable obscurer in 
NGC 5548 is a new breed of weakly-ionised, higher-velocity 
outflowing gas, which is much closer to the black hole and 
extends to the BLR. As reported in Kaastra et al. (2014) the 
outflowing obscurer is likely to originate from the accretion 
disk. Based on the high outflow velocity of the obscurer, its 
short-timescale absorption variability, and its covering fractions 
of the continuum and the BLR, the obscurer is in close proximity 
to the central source, and its geometry extends from near the 
disk to outside the BLR. 

In this work we present a broadband spectral analysis of the 
NGC 5548 data. The structure of the paper is as follows. Section 

2 gives an overview of our multi-satellite campaign. In Sect. 

3 we present lightcurves of NGC 5548 constructed at various 
energies from near-infrared (NIR) to hard X-rays. In Sect. 4 we 
explain the required steps in determining the spectral energy 
distribution (SED) of NGC 5548. In Sect. 5 we examine the soft 
X-ray excess in NGC 5548 and present an appropriate model 
for it. In Sec. 6 we describe the modelling of the broadband 
continuum in unobscured and obscured epochs and present our 
results. The thermal stability curves, corresponding to various 
ionising SEDs, are presented in Sect. 7. We discuss all our 
findings in Sect. 8 and give concluding remarks in Sect. 9. The 
processing of the data from all the instruments is described in 
Appendix A. 
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The spectral analysis and modelling, presented in this work, 
were done using the SPEX^ package (Kaastra et al. 1996) ver¬ 
sion 2.05.02. We also made use of tools in NASA’s HEASOFT^ 
v6.14 package. The spectra shown in this report are background- 
subtracted and are displayed in the observed frame, unless 
otherwise stated in the text. We use C-statistics (Cash 1979) 
for spectral fitting (unless otherwise stated) and give errors 
at Icr (68%) confidence level. The redshift of NGC 5548 is 
set to 0.017175 (de Vaucouleurs et al. 1991) as given in the 
NASA/IPAC Extragalactic Database (NED). The adopted cos¬ 
mological parameters for luminosity computations in our mod¬ 
elling are Hq - 70 km s“* Mpc^*, Qa = 0.70 and Q.,„ - 0.30. 

2. Multi-wavelength campaign on NGC 5548 

At the core of our campaign in summer 2013 (22 June to 1 
August), there were 12 XMM-Newton observations, of which 
five were taken simultaneously with HST COS, four with 
INTEGRAL and two with NuSTAR observations. Throughout 
our campaign and beyond. Swift monitored NGC 5548 on a 
daily basis. There were also optical monitorings with WO and 
OCA. The summer XMM-Newton observations were followed 
by three Chandra LETGS observations taken in the first half 
of September 2013, one of which was taken simultaneously 
with NuSTAR. The September observations were triggered upon 
observing a large jump in the X-ray flux from our Swift mon¬ 
itoring. However, due to scheduling constraints by the time the 
triggered Chandra observations were made, the week-long peak 
of high X-ray flux was just missed and the tail end of the flare 
was caught. Nonetheless, the X-ray flux was still higher than 
during the XMM-Newton observations and improved LETGS 
spectra were obtained. During this autumn period (Sep-Nov 
2013), NGC 5548 was not visible to XMM-Newton and thus no 
XMM-Newton observations could be triggered. 

A few months later, two more XMM-Newton observations 
were taken, one in December 2013 and the other in Eebruary 
2014. The former observation was simultaneous with HST COS 
and NuSTAR observations, and the latter close in time to an 
INTEGRAL observation. The timeline of all the observations 
in our campaign is displayed in Eig. 1 and the observation logs 
are provided in Table 1 . In Appendix A, we describe the obser¬ 
vations made by the instruments of each observatory and give 
details on their data reduction and processing. 

3. Lightcurves of NGC 5548 

Here we present lightcurves of NGC 5548, from NIR to hard 
X-rays, obtained from the eight observatories used during our 
summer 2013 campaign. 

3.1. X-ray lightcurves 

Eigure 2 shows the X-ray lightcurves of NGC 5548, from 
the Swift, XMM-Newton, NuSTAR, INTEGRAL and Chandra 
observations, with the average energy flux of each observation 
calculated over four X-ray energy bands from soft to hard 
X-rays. The X-ray fluxes have been calculated by fitting each 
dataset over the required energy band with the model described 
later on in Sect. 6. 


* http://www.sron.nl/spex 
^ http://heasarc.nasa.gov/lheasoft 


Eor the XMM-Newton observations, the lowest X-ray fluxes 
were observed in Obs. 1: Eo. 3^2 kev = 1-40 x 10“^^ erg cm“^ s“' 
and E 2 -iokeV = 1-51 x 10“" erg cm“^ s“'. These fluxes, 
compared to those from the unobscured 2000 and 2001 ob¬ 
servations, are smaller by factors ranging from 17-27 (for 
0.3-2 keV) and 2-3 (for 2-10 keV). Later on, in the summer 
XMM-Newton campaign, the X-ray fluxes increased from their 
minimum in Obs. 1 by a factor of 2.8 (0.3-2 keV) and 2.3 
(2-10 keV), peaking at XMM-Newton Obs. 4. The X-ray flux 
during the summer XMM-Newton campaign was low relative 
to other obscured epochs, as seen from the long-term Swift 
monitoring. The average fluxes of the summer 2013 XMM- 
Newton observations (Eo. 3 - 2 keV = 3.02 x 10“'^ ergcm^^ s“' 
and E 2 -iokeV = 2.77 x 10“'* erg cm“^ s“*) are smaller by a 
factor of 1.6 (for 0.3-2 keV) and 1.1 (for 2-10 keV) than the 
average fluxes for the whole obscured epoch observed with 
Swift between Eebruary 2012 and July 2014 (Mehdipour et al. 
in prep). 

After the end of the summer XMM-Newton observations 
(Obs. 1-12), Swift kept on monitoring NGC 5548 until the mid¬ 
dle of September 2013, when it passed out of Swift’s visibility 
window. However, prior to that a sudden increase in the X-ray 
flux was observed with Swift which started on 24 August 2013 
(56528 in MJD) and lasted about one week. Compared to the av¬ 
erage flux of the summer XMM-Newton observations, the flux in¬ 
creased by a factor of 3.4 (for 0.3-2 keV) and 1.9 (for 2-10 keV). 
This resulted in the triggering of the Chandra LETGS observa¬ 
tions. Later on, the two final XMM-Newton observations (Obs. 
13 and 14) were taken in December 2013 and Eebruary 2014 
(these are outside the range of Eig. 2 lightcurves). Compared to 
the average XMM-Newton flux during the summer 2013 obser¬ 
vations, the Obs. 13 flux was slightly lower by a factor of 1.2 
(for 0.3-2 keV) and 1.1 (for 2-10 keV), and for Obs. 14 the flux 
was higher by a factor of 1.5 for 0.3-2 keV and about the same 
for 2-10 keV. A comprehensive analysis of the XMM-Newton 
lightcurves and their properties is reported in Cappi et al. (in 
prep). 

3.2. NIR/optical/UV lightcurves 

In Eig. 3 we show the optical/UV lightcurves from Swift 
UVOT, XMM-Newton OM, WO and OCA, taken in the I, R, 
V, B, U, UVWl, UVM2 and UVW2 Alters. Eigure 4 shows 
the HST COS fluxes variability at six narrow UV energy 
bands (which are free of spectral features and represented by 
their central wavelengths) between 1160 and 1793 A for the 
summer 2013 observations. To calculate the flux in broadband 
photometric Alters, one requires convolution of the response of 
the instrument and the transmission of the Alter as a function 
of wavelength, with the shape of the optical/UV spectrum, to 
convert from instrumental units to energy flux at an effective 
wavelength. Eor the shape of the optical/UV spectrum, rather 
than using default spectral shapes commonly incorporated in 
the data reduction softwares (such as standard stars spectra, 
or a power-law), we used the optical/UV spectrum model that 
we derived for NGC 5548 from modelling the OM and UVOT 
Optical and UV grism spectra and the HST COS continuum 
data (see Sect. 4 and Eig. 5). This approach allows for a more 
accurate flux calculation in the photometric Alters, taking into 
account the presence of strong AGN emission lines with a 
realistic continuum model. Eor instance, for the six Swift UVOT 
Alters, the custom-fit spectral model results in 5 to 15% flux 
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Summer/autumn 2013 campaign on NGC 5548 

15-Jun-13 l-Jul-13 15-JUH3 l-Aug-13 15-Aug-13 1-Sep-13 15-Sep-13 

OCA 
WO 
Chandra 
INTEGRAL 
NuSTAR 
HST/COS 
XMM-Newton 
Swift 

20 40 60 80 

Days (MJD) - 56455 

Winter 2013/2014 campaign on NGC 5548 


16-Dec-13 1-Jan-14 16-Jan-14 1-Feb-14 
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Fig. 1. Timeline of our multi-wavelength campaign on NGC 5548. The thickness of each rectangular symbol on the time axis is 
indicative of the length of that observation. The days in which Swift, WO and OCA made observations are indicated by crosses. 

Table 1. Observation log of the NGC 5548 campaign. For the Swift, OCA and WO monitorings, the Obs. numbers correspond to 
days in which observations were taken. For Swift, only recent observations taken in 2013-2014 are reported, including all monitoring 
programs during this period. The Swift lengths in ks are the total length of the observations in each year. For HST COS the span of 
each observation in ks is given. 


Observatory 

Obs. 

ID 

Start time (UTC) 
yyyy-mm-dd hh:mm 

Length 

(ks) 

Observatory 

Obs. 

ID 

Stai't time (UTC) 
yyyy-mm-dd hh:mm 

Length 

(ks) 

XMM-Newton 

1 

0720110301 

2013-06-22 03:53 

50.5 

HST COS 

1 

lc7001 

2013-06-22 13:25 

13.0 

XMM-Newton 

2 

0720110401 

2013-06-29 23:50 

55.5 

HST COS 

2 

lc7002 

2013-07-12 02:23 

14.2 

XMM-Newton 

3 

0720110501 

2013-07-07 23:28 

50.9 

HST COS 

3 

lc7003 

2013-07-24 16:43 

8.9 

XMM-Newton 

4 

0720110601 

2013-07-11 23:11 

55.5 

HST COS 

4 

lc7004 

2013-07-30 15:15 

16.0 

XMM-Newton 

5 

0720110701 

2013-07-15 22:56 

55.5 

HST COS 

5 

lc7005 

2013-08-01 03:21 

12.3 

XMM-Newton 

6 

0720110801 

2013-07-19 22:40 

56.5 

HST COS 

6 

lc7006 

2013-12-20 02:49 

13.0 

XMM-Newton 

7 

0720110901 

2013-07-21 22:32 

55.5 

NuSTAR 

1 

60002044002 

2013-07-11 09:50 

51.6 

XMM-Newton 

8 

0720111001 

2013-07-23 22:24 

55.5 



60002044003 

2013-07-12 00:10 

52.2 

XMM-Newton 

9 

0720111101 

2013-07-25 22:15 

55.5 

NuSTAR 

2 

60002044005 

2013-07-23 14:25 

97.2 

XMM-Newton 

10 

0720111201 

2013-07-27 22:06 

55.5 

NuSTAR 

3 

60002044006 

2013-09-10 21:25 

97.5 

XMM-Newton 

11 

0720111301 

2013-07-29 21:58 

50.4 

NuSTAR 

4 

60002044008 

2013-12-20 08:30 

98.1 

XMM-Newton 

12 

0720111401 

2013-07-31 21:49 

55.5 

Chandra 

I 

16369 

2013-09-01 00:01 

29.7 

XMM-Newton 

13 

0720111501 

2013-12-20 14:01 

55.3 

Chandra 

2 

16368 

2013-09-02 10:33 

67.5 

XMM-Newton 

14 

0720111601 

2014-02-04 09:33 

55.5 

Chandra 

3 

16314 

2013-09-10 08:17 

122.0 

INTEGRAL 

1 

10700010001 

2013-06-29 21:34 

100.0 

Swift iiony 

1-160 

h 

2013-01-04 00:24 

326.6 

INTEGRAL 

2 

10700010002 

2013-07-11 21:13 

102.0 

Swift 

161-291 

d 

2014-01-02 14:53 

182.5 

INTEGRAL 

3 

10700010003 

2013-07-15 23:31 

106.5 

OCA" 

1-27 


2013-05-20 03:17 

e 

INTEGRAL 

4 

10700010004 

2013-07-23 03:54 

128.9 

WO (2013)f 

1-93 

- 

2013-06-02 07:22 

f 

INTEGRAL 

5 

11200110003 

2014-02-09 10:00 

38.2 

WO (2014) s 

94-150 

- 

2013-12-16 14:38 

g 


“ The monitoring in 2013 ended on 2013-12-31 18:19. 

* The target IDs of NGC 5548 in 2013: 30022, 80131, 91404, 91711, 91737, 91739, 91744, 91964. 

The Swift monitoring up to 2014-07-01 00:00 is reported here. The Swift monitoring of NGC 5548 is currently ongoing in 2014-2015. 
'' The Swift target IDs of NGC 5548 in 2014: 30022, 33204, 91964. 

The OCA monitoring ended on 2013-07-25 01:13. Observations taken in the B, V, R filters, with 150 s exposure in each filter. 

^ The WO monitoring ended on 2013-09-24 04:59. Observations taken in the B, V, R, I filters, with 300 s exposure in each filter. 

^ The WO monitoring ended on 2014-04-14 08:25. Observations taken in the B, V, R, I filters, with 300 s exposure in each filter. 
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Fig. 2. X-ray lightcurves of NGC 5548 from our summer 2013 campaign. The lightcurves show the observed flux and are displayed 
between 7 June and 14 Sep 2013. The XRT data are shown as blue circles, EPIC-pn as magenta stars, LETGS as dark yellow 
squares, NuSTAR as red diamonds and INTEGRAL data as green squares. The vertical dashed lines indicate the interval between 
XMM-Newton Obs. 1 and 12. 



20 40 60 80 

Days (MJD) - 56450 


20 40 60 80 

Days (MJD) - 56450 


Fig. 3. NIR, optical and UV lightcurves of NGC 5548 taken with various Alters from our summer 2013 campaign. The lightcurves 
show the observed energy flux and are displayed between 7 June and 14 Sep 2013. The UVOT data are shown as blue circles, 
OM as magenta stars, WO as red circles, and OCA data as green squares. The vertical dashed lines indicate the interval between 
XMM-Newton Obs. 1 and 12. 
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Fig. 4. HST COS fluxes of NGC 5548 at different UV wavelengths from summer 2013 observations. The lightcurves {left panel) 
and spectra {right panel) show the observed energy flux obtained from narrow energy bands which are free of spectral features. 
The displayed range of the lightcurves is between 16 June and 6 August 2013, with the vertical dashed lines indicating the interval 
between AMM-Aewfon Obs. 1 and 12. 


improvement in our knowledge of the continuum flux in each 
filter. The calculation of the COS fluxes at the narrow UV 
energy bands is described in Appendix A.6. In the lightcurves 
of Fig. 3, the UVOT and OM fluxes have been calculated at 
the following wavelengths for each filter: V (5402 A), B (4329 
A), U (3501 A), UVWl (2634 A), UVM2 (2231 A), UVW2 
(2030 A). For OCA and WO the fluxes have been given at the 
following wavelengths for each filter: I (8060 A), R (7000 A), V 
(5500 A), B (4330 A). 

As illustrated in the Fig. 3 lightcurves, the NIR/optical/UV 
fluxes were more or less continuously increasing throughout 
our summer XMM-Newton campaign (Obs. 1-12) and beyond, 
until the end of our Swift monitoring in middle of September 
2013, when the source was not visible to Swift for a couple of 
months. In fact the optical/UV fluxes recorded by Swift at mid 
September were the highest observed in all of the Swift monitor¬ 
ing of NGC 5548 spanning from 2005 to 2014, which also coin¬ 
cided with our triggered Chandra LETGS observations. Figure 4 
shows that the COS UV fluxes increased from their minimum at 
XMM-Newton Obs. 1 (COS Obs. 1) and peaked at XMM-Newton 
Obs. 8 (COS Obs. 3), followed by a decrease afterwards. The 
COS flux change from the first observation gets larger towards 
higher UV energies: factor of 1.8 at 1793 A to factor of 2.3 at 

1160 A. 

4. Determination of the continuum emission 

In order to establish the actual intrinsic continuum in the 
NIR/optical/UV part of the SED, the following absorption 
effects and emission contributions must be taken into account: 
(1) Galactic interstellar reddening; (2) Host galaxy stellar emis¬ 
sion; (3) Emission lines from the broad-line region (BLR) and 
narrow-line region (NLR); (4) Eeii blended emission feature; 
(5) Balmer continuum emission. Moreover, in the X-ray band 
the following must be taken into account: (6) Galactic X-ray 
absorption; (7) Absorption by the obscurer; (8) Absorption by 
the warm absorber; (9) Soft X-ray emission lines. 

In the following subsections, we describe the modelling of 
the above NIR/optical/UV components (1-5) in Sects. 4.1-4.5 


and the modelling of the X-ray components (6-9) in Sects. 
4.6-4.9. We have modelled the grism spectra from OM and 
UVOT together with simultaneous HST COS continuum mea¬ 
surements to establish the contribution of each component and 
consequently correct the data taken in the photometric filters. 
Eor fitting the grism spectra in SPEX we used minimisa¬ 
tion instead of C-statistics (used for X-ray spectra), as the op¬ 
tical/UV spectra have sufficient counts per bin. In Eig. 5 the 
best-fit model to the stacked OM and UVOT grism spectra and 
HST COS continua obtained from our simultaneous 2013 ob¬ 
servations are shown. In Eig. 6 we display all the stacked data 
from our summer 2013 campaign. The figure shows how much 
the aforementioned NIR/optical/UV corrections modify the data 
from observed (left panel) to corrected (right panel); the emerg¬ 
ing shape of the NIR/optical/UV SED after the corrections is 
remarkable and exhibits a thermal disk spectrum. In Eig. 6 the 
displayed X-ray data are only corrected for Galactic interstellar 
absorption and still include the effects of heavy absorption by the 
obscurer and the ‘de-ionised’ warm absorber (i.e. less ionised 
warm absorber due to the shielding of the ionising radiation by 
the obscurer). Similarly, in Eig. 7 the corrected optical/UV (OM) 
and X-ray (EPIC-pn and RGS) data from 2000, 2001 and av¬ 
erage summer 2013 are displayed, demonstrating that whilst in 
2013 the optical/UV continuum is higher than during the un¬ 
obscured epochs, the soft X-rays are heavily suppressed by the 
obscurer. 


4.1. Galactic interstellar reddening 


To correct the optical/UV fluxes for interstellar reddening in our 
Galaxy, the reddening curve of Cardelli et al. (1989) was used, 
including the update for near-UV given by O’Donnell (1994). 
The ebv model in SPEX applies this de-reddening to the data. Eor 
NGC 5548 the colour excess is E{B - V) - 0.02 mag (Schlegel 
et al. 1998). The scalar specifying the ratio of total to selective 
extinction Ry = AvlE{B - V) was set to 3.1. In Eig. 5, the red¬ 
dened continuum model is shown in dotted black line and the 
de-reddened one in dashed black line. 
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4.2. Host galaxy stellar emission 

To take into account the contribution of starlight from the host 
galaxy of NGC 5548 in our NIR/optical/UV data, we pro¬ 
duced an appropriate model for inclusion in our spectral mod¬ 
elling. Bentz et al. (2009,2013) have determined the host galaxy 
flux at an optical wavelength for a sample of AGN (includ¬ 
ing NGC 5548) using HST. However, for our observations of 
NGC 5548, this host galaxy flux was re-calculated to take into 
account the 5 arcsec radius circular aperture used in our pro¬ 
cessings. For the HST F550M medium-band V Alter, the host 
galaxy flux in a 5 arcsec radius circular aperture was determined 
to be 6.2 X 10“^^ erg cm“^ s“' with an uncertainty of ~ 10% 
(Misty Bentz, private communication). The effective wavelength 
of this flux measurement is at an observed wavelength of 5580 
A. Then, in order to calculate the host galaxy spectrum at the 
other wavelengths, we used a template model spectrum and nor¬ 
malised it to NGC 5548 host galaxy flux at 5580 A. A 5 arcsec 
aperture takes in only the innermost few kpc of the host galaxy, 
and so the galaxy bulge template of Kinney et al. (1996) was 
adopted. In Fig. 5, the contribution of the host galaxy stellar 
model is displayed as the solid green line. 


4.3. Emission lines from BLR and NLR 

To take into account emission lines produced from the BLR and 
NLR, we modelled the OM Optical grism and UVOT UV grism 
spectra, covering a range from 1895 to 6850 A. The simultane¬ 
ity of the grism spectra with COS is important as it enables 
us to model the grism spectra and COS continuum points to¬ 
gether. This wider energy band in the optical/UV helps in es¬ 
tablishing the underlying optical/UV continuum and modelling 
the BLR and NLR lines. For the underlying optical/UV contin¬ 
uum, the Comptonisation component comt in SPEX (explained 
later in Sects. 5 and 6) was used. The broad and narrow emis¬ 
sion lines were modelled using Gaussian line profiles. For broad 
lines, multiple Gaussian components with different widths, but 
the same wavelength were added until a good fit was obtained 
for each line. The best fit obtained for the BLR and NLR lines 
was saved as a spectral model component which was included 
in all our broadband modelling later in Sect. 6, to take into ac¬ 
count the contribution of these lines in the photometric Alters. 
The flux contribution of all the BLR/NLR lines in the grism en¬ 
ergy range shown in Fig. 5 (solid black lines) was found to be 
about 7.94 x 10“'^ erg cm“^ s“'. 


4.4. Fe II blended emission feature 

The Fen in the BLR produces several thousands of transi¬ 
tions, which result in a blended and complex spectrum between 
-2000^000 A (see e.g. Netzer & Wills 1983; Wills et al. 1985). 
In order to take into account the Fe ii contribution in our opti¬ 
cal data we took the model calculated by Wills et al. (1985), 
which is convolved to an intrinsic width of 2500 km s '. This 
was then imported into SPEX as a spectral model component and 
convolved with the resolution of the grisms. The normalisation 
scaling factor of the component was left as a free parameter to 
allow for fitting the total Fe ii flux in NGC 5548. The flux contri¬ 
bution by Fe ii in NGC 5548 was found to be about 2.79 x 10“'^ 
erg cm “ s '. This component is shown in Fig. 5 as part of the 
magenta coloured feature. 


4.5. Balmer continuum feature 

Apart from Fe ii there is another significant emission component 
in the ~2000^000 A region of AGN spectra: the Balmer 
continuum (see e.g. Grandi 1982). The Balmer continuum for 
an optically-thin emission with constant electron temperature 
is given by F®*- = ^BEg-/i(y-yBE)AA (Qj-^ndi 1982), where F®'' 
is the Balmer continuum flux at frequency v, F®® is the flux at 
the Balmer edge vbe, h the Planck constant, k the Boltzmann 
constant and the electron temperature. The Balmer edge is at 
theoretical wavelength of 3646 A. We used this model to take 
into account the contribution of the Balmer continuum in our 
optical data. The model was convolved with an intrinsic Doppler 
velocity broadening cr„ and the resolution of the grisms in SPEX. 
We fitted this model to the grism data for various combinations 
of Fe and cr„ values and found the best-fit Balmer continuum 
component has Fe ~ 8000 K and cr„ ~ 10000 km s '. We note 
that the purpose of fitting the Balmer continuum in this work is 
to obtain its total flux contribution to the optical data in order 
to establish the underlying continuum component, so we do 
not delve into the Fe and cr„ parameters here which carry large 
uncertainties. The flux contribution by the Balmer continuum 
was found to be about 5.75 x 10“'^ erg cm“^ s“'. 

The Balmer continuum, together with Fen emission, form a 
broad and blended feature in AGN spectra dubbed the ‘small- 
blue-bump’. This is displayed in Fig. 5, in which our best-fit 
model to the grism spectra and HST COS continua points are 
shown. The small-blue-bump is the feature in magenta colour 
between ~2000^000 A; beneath that the Balmer continuum 
model (without Fe ii emission) is displayed in dark yellow. 


4.6. Galactic Interstellar X-ray absorption 

The effects of the Galactic neutral absorption in the interstel¬ 
lar medium are included in our modelling by applying the hot 
model in SPEX. This model calculates the transmission of gas in 
collisional ionisation equilibrium. For a given temperature and 
set of abundances, the model calculates the ionisation balance 
and then determines all the ionic column densities by scaling 
to the prescribed total hydrogen column density Na- The trans¬ 
mission includes both continuum and line opacity. To mimic the 
transmission of a neutral plasma in collisional ionisation equilib¬ 
rium (such as the interstellar medium of our Galaxy), the temper¬ 
ature of the plasma is set to 0.5 eV. In our modelling the Galactic 
HI column density in the line of sight to NGC 5548 was fixed to 
Vh = 1.45 X 10^® cm“^ (Wakker et al. 2011) with Lodders et al. 
(2009) abundances. 


4.7. Absorption by the obscurer 

To take into account the newly-discovered X-ray obscuration in 
NGC 5548 we adopted the same model found by Kaastra et al. 
(2014). In this model the obscurer consists of two ionisation 
phases, each modelled with an xabs component in SPEX. The 
xabs model calculates the transmission through a slab of pho- 
toionised gas where all ionic column densities are linked in a 
physically consistent fashion through the CLOUDY photoionisa¬ 
tion model. The fitted parameters of an xabs model component 
are the ionisation parameter (^), the equivalent hydrogen column 
density (Vh), the covering fraction C/ of the absorber, its flow 
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Fig. 5. De-reddened optical/UV spectrum of NGC 5548. The displayed data are from simultaneous UVOT UV grism spectra (shown 
in light blue) at 1896-3414 A, OM Optical grism spectra (shown in dark blue) at 3410-6836 A and the six HST COS continuum 
points (shown as purple squares) between 1160-1793 A. The displayed data are the average of the contemporaneous OM, UVOT 
and HST COS observations taken in summer 2013 and are corrected for Galactic reddening using the model described in Sect. 4.1. 
The best-fit model (described in Sect. 4) is shown in red and various components contributing to the model are also displayed. The 
underlying continuum model (comt) is the dashed black curve. The dotted black curve is the reddened version of the continuum 
to illustrate the correction for reddening. The broad and narrow emission line components (Sect. 4.3) are shown at the bottom and 
some of the prominent emission lines are labelled. The broad feature in magenta is the ‘small-blue-bump’: blended Feii emission 
(Sect. 4.4) with Balmer continuum (Sect. 4.5). The model for the Balmer continuum alone is shown in dark yellow below the 
small-blue-bump. The contribution from the host galaxy of NGC 5548 (Sect. 4.2) is shown in green. 




Fig. 6. Average multi-wavelength data of NGC 5548 from our summer 2013 campaign. The data are displayed before {left panel) 
and after {right panel) corrections for Galactic reddening (Sect 4.1), host galaxy stellar emission (Sect 4.2), emission lines from the 
BLR and NLR (Sect 4.3), blended Fe ii emission (Sect. 4.4), Balmer continuum (Sect. 4.5) and Galactic interstellar X-ray absorption 
(Sect 4.6). The spectra have been binned for clarity of presentation. The displayed X-ray data include the effects of heavy absorption 
by the obscurer and the warm absorber. The shape of the NIR/optical/UV data after the corrections {right panel) is remarkable and 
exhibits a thermal disk spectrum. 
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Fig. 7. Optical/UV (OM) and X-ray (RGS and EPIC-pn) data of 
NGC 5548 from unobscured (2000, 2001) and obscured (2013) 
epochs. The optical/UV data have been corrected for the effects 
described in Sects. 4.1-4.5. The X-ray spectra have been cor¬ 
rected for Galactic absorption (Sect. 4.6), but not for absorption 
by the obscurer and the warm absorber. The displayed X-ray data 
between 0.3-1.4 are RGS and at higher energies EPIC-pn. The 
X-ray spectra have been binned for clarity of presentation. In 
the 2013 data, the soft X-ray spectrum is the lowest one and the 
optical/UV data the highest one. 


V and turbulent cr^ velocities. The ionisation parameter ^ (Tarter 
et al. 1969) is defined as 



where L is the luminosity of the ionising source over the 1- 
1000 Ryd (13.6 eV to 13.6 keV) band in erg s ', tin the hy¬ 
drogen density in cm“^ and r the distance between the ionised 
gas and the ionising source in cm. The first component of the 
obscurer covers about 86% of the central X-ray emitting region, 
with logf = -1.2 and Nh = 1.2 x 10^^ cm“^. The second com¬ 
ponent of the obscurer covers 30% of the X-ray source and is 
almost neutral (log^ = -4.0) with Nn - 9.6 x 10^^ cm“^. The 
parameters of the obscurer were fixed to those of Kaastra et al. 
(2014) in our broadband spectral modelling of the stacked sum¬ 
mer 2013 data. Of course for the unobscured epochs (2000 and 
2001 XMM-Newton data), the obscurer components were ex¬ 
cluded from our model by setting C / of both obscurer compo¬ 
nents to zero. 

4.8. Absorption by the warm absorber outflows 

Since the obscurer is located between the central ionising source 
and the warm absorber, it prevents some of the ionising radia¬ 
tion from reaching the warm absorber. Thus the different phases 
of the warm absorber become less ionised, resulting in more X- 
ray absorption than when NGC 5548 was unobscured. In our 
modelling of the obscured data we used the de-ionised warm ab¬ 
sorber model obtained by Kaastra et al. (2014), which consists 
of six different phases of photoionisation. Each phase is repre¬ 
sented by an xabs component in SPEX with Nh, flow v and 
turbulent cr„ velocities kept frozen throughout our modelling to 
the values given in Kaastra et al. (2014). The warm absorber 
phases have Nh ranging from 2 to 57 xlO^" cm“^, log^ from 


0.33 to 2.67, (T„ from 20 to 210 km s“' and outflow velocities 
ranging from about 250 to 1220 km s“'. Eor the 2000 and 2001 
XMM-Newton data, when the warm absorber would have been 
exposed to normal ionising flux, parameters of the warm ab¬ 
sorber were taken from a re-analysis of their RGS data, using 
the same method as in Kaastra et al. (2014). This work will be 
reported in a forthcoming paper on our campaign (Ebrero et al. 
in prep), where the long-term variability of the warm absorber is 
going to be presented. 

4.9. Soft X-ray emission lines 

The 2013 RGS data clearly show the presence of several narrow 
emission lines and radiative recombination continua from pho- 
toionised gas (Kaastra et al. 2014). These emission features are 
more apparent in 2013 as the soft X-ray continuum is suppressed 
due to obscuration, but they are also present at earlier epochs. In 
order to take into account their flux contribution in our contin¬ 
uum modelling we used the same model and parameters reported 
in Kaastra et al. (2014). The contribution of these lines to the to¬ 
tal 0.3-2.0 keV flux is 1.4% (in 2000), 0.8-1.3% (in 2001) and 
7.7% (in 2013). So whilst in the unobscured epochs their contri¬ 
bution is rather small, in obscured epochs they contribute more 
to the soft X-ray flux. In another forthcoming paper on our cam¬ 
paign by Whewell et al. (in prep), a detailed study of the soft 
X-ray emission features as detected by RGS in NGC 5548 is 
reported. 

5. The soft X-ray excess in NGC 5548 

The AGN ‘soft X-ray excess’ is an excess continuum emission 
above the intrinsic X-ray power-law at the soft X-ray energies 
(below ~ 2 keV). Since its discovery by Singh et al. (1985) in 
HEAO-1 observations of Mrk 509 and also by Amaud et al. 
(1985) in EXOSAT observations of Mrk 841, it has been ob¬ 
served in many Seyfert-1 AGN. Eor example, in the Catalogue 
of AGN In the XMM-Newton Archive (CAIXA, Bianchi et al. 
2009), this component was commonly found in about 80% of 
the AGN. In NGC 5548 the soft excess has been previously 
detected (e.g. Kaastra & Barr 1989; Magdziarz et al. 1998) and 
its presence in the unobscured state is apparent. Eigure 8 shows 
the EPIC-pn spectrum obtained in 2000. The continuum above 
3.0 keV has been fitted with a Galactic-absorbed power-law 
component (T ~ 1.8), including a reflection component for 
modelling the Ee Ka line. The best-fit is then extrapolated 
to lower energies, displaying the presence of the soft excess 
below about 2 keV. However, the soft excess in NGC 5548 only 
stands out in data before the X-ray obscuration began (i.e. < 
2011). Because of the newly discovered obscurer (Kaastra et al. 
2014), the soft X-ray flux is heavily absorbed in 2013. This 
demonstrates that the soft excess emission is produced in a small 
region interior to the obscurer (e.g. the BLR) and is not from a 
large-scale scatterer. Yet, the soft X-ray band shows significant 
variability which can be attributed to both the obscurer and/or 
soft excess component. Therefore, for any variability study of 
the obscurer, the variability of the soft excess component also 
needs to be established. This means an appropriate model for 
the soft excess is required which we describe below. 

There are different interpretations and models in the litera¬ 
ture for the nature of the enigmatic soft excess emission in AGN 
(see Sect. 8.1). Eor the case of NGC 5548, evidence (described 
below) suggests that the soft X-ray excess is primarily part 
of an optical/UV/soft X-ray continuum component, produced 
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Fig. 8. EPIC-pn spectrum of NGC 5548 from 2000, fitted above 
3.0 keV with a Galactic-absorbed power-law (E ~ 1.8), includ¬ 
ing a reflection component for the Ee Ka. The ht is extrapolated 
to lower-energies displaying the presence of a soft excess be¬ 
low about 2 keV. The ht residuals, (observed-model)/model, are 
displayed in the bottom panel. 


by ‘warm Comptonisation’: up-scattering of the seed disk 
photons in a Comptonising corona, which is distinct from the 
one responsible for the hard X-ray power-law (i.e. the hot 
corona), and is characterised by a lower temperature and higher 
optical depth (e.g. Magdziarz et al. 1998). Eirstly, without any 
modelling, the unobscured NGC 5548 data (Eig. 9, top panel) 
show that as the UV Hux increases, the ratio of soft to hard 
X-ray hux (‘softness’ ratio) also increases. This indicates the 
UV and soft excess are related to each other. In Eig. 9 (top 
panel) the X-ray softness ratio is plotted versus the UV hux for 
the unobscured data (XMM-Newton 2000 and 2001 and Swift 
2007), and also for the obscured XMM-Newton data from the 
average of the summer 2013 observations. The softness ratio 
corresponds to the observed Hux, without any modelling and 
corrections for absorption. Remarkably, for the obscured epoch 
in 2013, the softness ratio is at its lowest level even though the 
UV hux is at its highest. The reason for this sharp drop in the 
softness ratio is the presence of the heavy soft X-ray absorption 
by the obscurer in 2013 (see Eig. 7). 

The link between the UV and soft excess hux was also 
seen in Mrk 509 (Mehdipour et al. 2011). Erom broadband 
modelling of its continuum, Mehdipour et al. (2011) and 
Petrucci et al. (2013) argued warm Comptonisation to be the 
most reasonable explanation for the soft excess. Eurthermore, 
warm Comptonisation has already been suggested for the soft 
excess in NGC 5548, in one of the hrst papers investigating this 
interpretation, by Magdziarz et al. (1998). Erom broad-band 
spectral analysis of NGC 5548, using data from lUE, Ginga, 
ROSAT, CGRO and OSSE, Magdziarz et al. (1998) found that 
the soft excess requires a separate continuum component which 
can be htted by Comptonisation of thermal photons from a cold 
disk (Tseed ~ 3.2 eV), in a warm (~ 0.3 keV), optically thick 
(t ~ 30) corona. The authors also found the optical/UV and soft 
X-ray Huxes in NGC 5548 to be closely correlated. Therefore 
considering the above evidence we adopt warm Comptonisation 
as the appropriate model for the soft excess. 



Fig. 9. Top panel: Ratio of observed soft (0.3-2.0 keV) to hard 
(2.0-10 keV) X-ray hux plotted versus the observed UVW2 
(2030 A) Hux. Bottom panel: Intrinsic luminosity of the soft 
excess component between 0.3-2.0 keV plotted versus the ob¬ 
served UVW2 (2030 A) Hux, and htted with the power-law func¬ 
tion given in Eq. 2. 

In Eig. 9 (bottom panel) we show the soft excess luminosity 
in the 0.3-2.0keV band versus the UV Hux, showing the relation 
between the soft excess and UV. The soft excess luminosity is 
calculated using the model described in Sect. 6 (i.e. corrected 
for obscurer, warm absorber and Galactic absorption). The data 
in Eig. 9 (bottom panel) have been htted with a function which 
is found to be a power-law, given by 

Tsoft excess = 1.261 (Eyvwz)^'*^^ (2) 

where Lsoft excess is the intrinsic luminosity of the soft excess 
component in the 0.3-2.0 keV band in units of 10^^ erg s * and 
Euvw 2 is the UVW2 hlter Hux in units of erg cm“^ s“* A“'. This 
Hux in the UVW2 hlter is at the wavelength of 2030 A for both 
UVOT and OM and is the observed hux without any modihca- 
tions such as de-reddening, but takes into account the optical/UV 
spectrum shape of NGC 5548 (shown in Eig. 5) in the UVW2 hl¬ 
ter bandpass. We further discuss the soft excess in Sect. 8.1. 

6. Broadband continuum modelling of NGC 5548 

Here we describe our continuum modelling of the NGC 5548 
data. The archival XMM-Newton data from 2000 and 2001 
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and the stacked 2013 data from our campaign have been used. 
The archival data are from epochs before the source became 
obscured and its spectrum more complex, and so are useful 
in establishing the underlying continuum of NGC 5548. The 
stacked 2013 data refer to stacked XMM-Newton data (Ohs. 1- 
12) and the corresponding simultaneous NuSTAR, INTEGRAL, 
HST COS, WO (I and R filters) observations collected during 
the summer 2013 campaign. The stacked observations are all 
from a period when the X-ray flux was at a stable low level 
with only moderate variability between observations, such that 
it does not have an adverse effect on the results obtained from 
the stacked spectrum. The stacking of spectra was performed 
using the mathpha, addrmf, and addarf tools of the HEASOFT 
package, except for RGS where the stacking was done using 
the procedure described in Appendix A. 1 . All the X-ray spectra 
used in our fitting have been optimally binned in SPEX (the obin 
command) taking into account the instrumental resolution and 
count statistics of the source as described in the SPEX manual. 
In our broadband spectral modelling we included all the model 
components discussed in Sect. 4, which take into account: 
Galactic reddening, host galaxy stellar emission, BLR and NLR 
emission lines, blended Lett and Balmer continuum. Galactic 
X-ray absorption, the warm absorber, soft X-ray emission lines 
and Anally absorption by the obscurer (for the 2013 data). 

As explained in Sect. 5, warm Comptonisation is a reason¬ 
able and likely model for the soft excess in NGC 5548. So in our 
broadband modelling of the continuum, we applied the comt 
model in SPEX to simultaneously fit the NIR/optical/UV and soft 
X-ray data. The comt model is based on the Comptonisation 
model of Titarchuk (1994), but with improved approximations 
for the parameter /3{t) which characterises the photon distri¬ 
bution over the number of scatterings which the soft photons 
undergo before escaping the plasma (for more details on this 
see the SPEX manual). The up-scattering Comptonising plasma 
was set to have a disk geometry and its parameters are: the 
temperature of the seed photons Tseed, the electron temperature 
of the plasma Tg, its optical depth t, and normalisation of the 
component. We note the parameterisation of the warm corona 
is model-dependent and different Comptonisation models can 
lead to slightly different T and t values. All models however 
give qualitatively the same results, i.e. low-T and high-r values 
and reproduce well the characteristic spectral shape of the soft 
excess. 

In addition to the comt component which models the 
NIR/optical/UV continuum and the soft X-ray excess, a 
power-law component (pow in SPEX) was also included in 
our model. This component models the primary hard X-ray 
continuum which has a power-law spectral shape and is 
thought to be produced by inverse Compton scattering in an 
optically thin (t ~ 1), hot (Tg ~ 100 keV) corona (e.g. Sunyaev 
& Titarchuk 1980; Haardt & Maraschi 1993). The cut-off 
power-law is known to provide a reasonable approximation for 
the Comptonisation spectrum (e.g. Petrucci et al. 2000). The 
high-energy exponential cut-off of the power-law was set to 
400 keV in our modelling; for details of the high-energy cut-off 
parameter search in NGC 5548 see Ursini et al. (submitted). 
We note the seeming decline in about 50-80 keV in Lig. 10 
is background dominated and not statistically significant. The 
actual high-energy cut-off is at higher energies (~300^00 keV). 
At low energies, the power-law was smoothly broken before 
overshooting the energies of the seed photons from the disk. 
The break energy was fixed to the peak energy of the comt 


component. The precise energy, which is rather uncertain, 
was found to have negligible effect in the NIR/optical/UV 
modelling of the NGC 5548 data, as the comt component is 
overly dominant at these energies (see Lig. 10). 

To model the Le Ka line and the reflection spectrum at hard 
X-ray energies we used the reflection model refl in SPEX. This 
model, provided by Piotr Zycki, computes the reflected contin¬ 
uum and the corresponding Le Ka line from a constant density 
X-ray illuminated atmosphere. It incorporates the Compton- 
reflected continuum of Magdziarz & Zdziarski (1995) and the 
Le Kg' line model of Zycki & Czerny (1994). The parameters of 
the model were fixed to those reported in Kaastra et al. (2014) 
with no relativistic broadening. We fitted the normalisation and 
photon index T of its illuminating power-law in our modelling 
to determine the cold reflection from distant material. However, 
for the 2000 and 2001 observations (in which no hard X-ray 
NuSTAR and INTEGRAL data were available), L of refl was 
frozen to the value found from 2013 data (L = 2.2 + 0.1) and 
only the normalisation was fitted in order to make a good fit 
to the Le Ka. The power-law normalisations of the reflection 
component were found to be 4.4 + 0.4 (for 2000), 3.5 + 0.3 (for 
9 July 2001), 4.4 + 0.4 (for 12 July 2001) and 3.0 + 0.4 (for 
stacked 2013) in units of 10^^ photons s ' keV ' at 1 keV. 

In our broadband spectral modelling, Tseed, Tg, t and nor¬ 
malisation of comt, as well as T and normalisation of pow were 
fitted. At the last step in fitting the 2013 data, we tested if free¬ 
ing parameters of the obscurer from their initial values (Kaastra 
et al. 2014) provides any improvement to the fit. However, the 
obscurer parameters remained within 1-cr errors and the good¬ 
ness of the fit was not significantly improved. Therefore, for 
consistency we have kept the obscurer parameters fixed to the 
values reported by Kaastra et al. (2014). Ligure 10 shows our 
broadband model fit to the stacked 2013 data. The contributions 
from individual continuum components (comt, pow and refl) 
are also displayed in the figure. In Table 2 the best-fit parame¬ 
ters of our continuum model are given. The observed fluxes and 
intrinsic luminosities in various energy bands, calculated from 
our model, are given in Table 3. In Lig. 11 we compare the in¬ 
trinsic broadband continuum SED at unobscured and obscured 
epochs. In Lig. 12 we show the broadband SED in 2013 before 
and after absorption by the obscurer. 

7. Thermal stability analysis of photoionised gas 

Having developed a broadband continuum SED model for the 
unobscured (2000 and 2001) and obscured (summer 2013) 
epochs, we investigated the effects of different SEDs and X-ray 
obscuration on the ionisation balance of gas. This is important 
for investigating the impact of the ionising SED on the obscurer 
and subsequently, after heavy absorption, on the warm absorber 
outflows. 

The ionising SED determines the ionisation balance and 
thermal stability of photoionised gas, such as warm absorbers 
in AGN. X-ray photoionised plasma can be thermally unstable 
in certain regions of the ionisation parameter space. This can be 
investigated by means of producing the thermal stability curves 
(also referred to as S-curves or cooling curves). These are plots 
of either the temperature of the plasma, or the ionisation param¬ 
eter versus the pressure form of the ionisation parameter, E, 
introduced by Krolik et al. (1981). The parameter S, which de¬ 
scribes the ionisation equilibrium, is defined as E s F/iin c kT, 
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0.001 


0.1 1 
Observed Energy (keV) 


Fig. 10. Best-fit broadband continuum model for the stacked summer 2013 data of NGC 5548 as described in Sect. 6. All the 
displayed data, labelled in the figure, have been corrected for the emission and absorption effects described in Sect. 4, including 
obscuration and warm absorption, required for determining the underlying continuum. The soft X-ray emission features seen in 
the RGS are shown on the continuum model for display purposes. The comt warm Comptonisation component (modelling the soft 
excess) is shown in black dashed-line, with the power-law (pow) and reflection (refl) components in black dotted-line form. The 
broadband model spectrum is represented in a solid green line. Residuals of the fit, defined as (data - model) / model, are displayed 
in the bottom panel. 


Table 2. Best-fit parameters of the broadband continuum model for 2000, 2001 and stacked 2013 data of NGC 5548. 


Observation 

Warm Comptonisation (comt): 

T’ b -c 

^ seed e ‘ 

Norm'* 

Power-law (pow): 

T" Norm^ 

C-stat / Expected C-stat 

2000 (24 Dec) 

0.74 ± 0.03 

0.15 ±0.01 

22.7 ± 0.9 

6.3 ±0.5 

1.82 ±0.01 

5.94 ± 0.08 

1544/ 1353 

2001 (09 July) 

0.71 ±0.03 

0.14 ±0.01 

22.1 ±0.8 

5.2 ±0.4 

1.81 ±0.01 

7.28 ± 0.04 

1860/ 1372 

2001 (12 July) 

0.58 ± 0.05 

0.13 ±0.01 

24.3 ± 1.0 

6.4 ± 0.7 

1.88 ±0.01 

10.3 ± 0.08 

1790/ 1362 

2013 (Summer) 

0.80 ± 0.02 

0.17 ±0.01 

21.1 ±0.5 

6.0 ± 0.2 

1.61 ±0.01 

4.72 ±0.10 

2666 / 1664 


“ Temperature of the seed photons in eV. 
* Warm corona temperature in keV. 
Optical depth. 

10^^ photons s“'keV“*. 

Photon index. 

^ 10^' photons s“* keV“* at 1 keV. 
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Table 3. Observed fluxes and intrinsic luminosities of NGC 5548 for 2000, 2001 and stacked 2013 data. The fluxes (F) are given 
in units of 10 " erg cm“^ s ' and the intrinsic luminosities (L) are in 10"'^ erg s *. The observed fluxes include all the absorptions 
described in Sects. 4.6, 4.7 and 4.8. The host galaxy stellar emission (Sect. 4.2) is excluded from the reported flux and luminosities. 


Observation 

Eo ,3 2 keV 

ko.3-2 keV 

E2-IO keV 

7.2-10 kcV 

Eo. 001-10 keV 

7 . 0 , 001-10 keV 

7^1-1000 Ryd 

7.1-1000 Ryd 

2000 (24 Dec) 

2.34 

2.53 

3.14 

2.17 

13.1 

18.4 

6.50 

12.2 

2001 (09 July) 

2.72 

2.47 

3.87 

2.67 

13.6 

15.8 

8.04 

10.9 

2001 (12 July) 

3.74 

3.44 

4.90 

3.39 

16.5 

19.6 

10.5 

14.4 

2013 (Summer) 

0.30 

2.21 

2.80 

2.35 

10.5 

17.9 

4.09 

11.7 



Fig. 11. The broadband continuum SED models for the unob¬ 
scured (2000, 2001) and obscured (2013) epochs of NGC 5548, 
obtained in Sect. 6. 



u. I I 

Rest Energy (keV) 



Fig. 12. The SED model for the 2013 epoch, displayed be¬ 
fore (dark blue) and after (light blue) absorption of the con¬ 
tinuum by the obscurer. The ionising luminosity Li_iooo Ryd is 
11.7 X 10''^ erg s“* before obscuration and 3.5 x 10^^ erg s“* 
after obscuration. 


where F is the flux of the ionising source between 1-1000 Ryd 
(in erg cm“^ s '), k the Boltzmann constant, T the temperature 
and «H is the hydrogen density in cm“^. Here, Fjc is the pres¬ 
sure of the ionising radiation. Taking F - LjAnr^ and using 
^ s L/uh S can be expressed as 

L £ £ 

H = -;- = —2-~ 19222 - (3) 

Ajir~ «H ckT A-nckT T 


Fig. 13. Thermal stability curves for photoionised gas in 
NGC 5548, showing ^ plotted versus H (top panel) and temper¬ 
ature kT plotted versus S (bottom panel). They have been calcu¬ 
lated for SEDs (displayed in Eigs. 11 and 12) derived from our 
broadband modelling of the continuum as described in Sect. 6. 
The circles on the stability curves indicate the positions at which 
ionic column densities of Ee peak on the (H, f) and (H, T) phase 
space. To guide the eye, the location of Eeix, Eexvii, Eexxv 
have been indicated with open circles around filled circles. 


On the curve itself, the heating rate is equal to the cooling 
rate, so the gas is in thermal balance. To the left of the curve, 
cooling dominates over heating, whilst to the right of the curve, 
heating dominates over cooling. On branches of the S-curve 
which have positive gradient, the photoionised gas is thermally 
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stable. This means small perturbations upward in temperature 
increase the cooling, whereas small perturbations downward in 
temperature increase the heating. However, on the branches with 
negative gradient, the photoionised gas is thermally unstable, in 
this case {d{C — H)ldT)^ < 0, where C represents the complete 
set of cooling processes and H represents the complete set of 
heating processes. Therefore on the unstable branches, a small 
perturbation upward in temperature increases the heating rela¬ 
tive to the cooling, causing further temperature rise, whereas a 
small perturbation downward increases the cooling. The thermal 
stability curve is highly sensitive to the SED of the ionising 
source and to the elemental abundances of the photoionised gas. 

We have computed the cooling curves (Fig. 13) correspond¬ 
ing to the different ionising SEDs of NGC 5548 (Figs. 11 and 
12). To obtain S(7’) and S(^), we ran CLOUDY^ (Ferland et al. 
1998) version C13.01 for a photoionised gas with Fodders et al. 
(2009) abundances irradiated with the SEDs obtained in Sect 6. 
In Fig. 13 the points on the cooling curves where the Fe ions 
peak have been indicated. Interestingly, the highly-ionised ions 
of Fe fall on an unstable branch of the curves for all hve cases. 
Moreover, the shielding by the obscurer in NGC 5548 has a sig- 
nihcant impact on the ionisation balance and thermal stability of 
photoionised gas located further out and produces a large unsta¬ 
ble branch on the curves. We further discuss the cooling curves 
results in Sect. 8.3. 

8. Discussion 

8.1. The nature of the NIR/optical/UV continuum and the soft 
X-ray excess 

The NIR/optical/UV continuum of NGC 5548, modelled using 
simultaneous data from about 8060 A to 1160 A, can be mod¬ 
elled as being one single continuum component. This result has 
been found in both local modelling of the optical/UV grisms and 
HST COS data (Sect. 4 and Fig. 5), and broadband modelling 
of the continuum (Sect. 6 and Fig. 10). This single continuum 
component is modelled by a warm Comptonisation component 
(comt) which, as well as fitting the NIR/optical/UV continuum, 
also hts the soft X-ray excess emission. There are no additional 
continuum components detected in the NIR/optical/UV ener¬ 
gies. We emphasise the importance of fully taking into account 
various non-continuum processes (described in Sect. 4) which 
make a signihcant flux contribution in broadband photometric 
hlters. For example, emission from the blended Fe ii and Balmer 
continuum (Sects. 4.4 and 4.5) contributes to 27% of the U filter 
flux, and emission from the host galaxy (described Sect. 4.2) 
contributes to about 48% of the V hlter flux. The host galaxy 
contribution relative to the continuum becomes greater in the R 
and I photometric Alters, contributing to 70% and 79% of the 
observed flux, respectively. So if only photometric Alter data 
are fitted without taking into account the actual spectrum and 
secondary emission processes, one may erroneously attribute an 
‘excess emission’ to an additional continuum component. 

Warm Comptonisation has been found to be a viable expla¬ 
nation for the soft excess in some AGN based on broadband 
spectral modelling or timing analyses - see e.g. Magdziarz 
et al. (1998); Mehdipour et al. (2011); Petrucci et al. (2013); 
Jin et al. (2013); Noda et al. (2013); Di Gesu et al. (2014); 
Arevalo & Markowitz (2014); Matt et al. (2014). The soft 

^ http://www.nublado.org 


excess has also been attributed to reprocessed emission caused 
by relativistically-blurred, ionised-disk reflection in some AGN 
(in particular NFSyls) based on observations of relativistically 
broadened Fe profiles and time lags between soft and hard 
X-rays - see e.g. Ross & Fabian (1993, 2005); Crummy et al. 
(2006); Fabian et al. (2009, 2013); Walton et al. (2013). 
In this work, the correlation between UV and soft X-rays 
long-timescale variability and the results of our broadband 
continuum modelling of NGC 5548 data suggest the origin of 
the soft excess is consistent with the warm Comptonisation 
interpretation. This was made possible by the remarkable 
simultaneous data from our summer 2013 campaign covering 
a wide energy range from NIR (8060 A) to hard X-rays (200 
keV) and also the archival (2000, 2001 and 2007) data when the 
source was not obscured. Our results suggest the soft excess is 
primarily the tail of the Comptonised disk component extending 
from NIR to soft X-rays. We note that although the long-term 
variability of the soft excess over several years can be explained 
by changes in the strength of the Comptonised disk emission, on 
short timescales reprocessed emission from reflection can also 
possibly contribute to any rapid variability of the soft excess. 
Indeed, short-timescale lag between soft and hard X-rays has 
been detected in NGC 5548 while the source was unobscured 
(De Marco et al. 2013), which may be interpreted as some 
contribution to soft excess rapid variability from reprocessing. 
However, the combination of the hard X-ray data of NuSTAR 
and INTEGRAF with EPIC-pn spectroscopy of the Fe Ka line 
show the reflection spectrum and the Fe Ko- line energy and 
its narrow width are consistent with reflection from distant 
neutral material (Ursini et al. submitted; Cappi et al. in prep). 
From our best-fit model the contribution of the neutral refl 
component to the soft excess is only ~ 1%. Although there is no 
spectroscopic evidence of an ionised reflection in NGC 5548, 
by invoking an additional ionised refl component to our model 
with log^ = 2.0 and relativistic blurring, the contribution of the 
ionised reflection component to the soft excess is estimated to 
be about 2-4%, with the rest of the soft excess emission being 
produced by the warm Comptonisation. 

The direct relation between the UV flux and the strength of 
the soft excess (Fig. 9) suggests they are related and primarily 
part of the same emission component driving the UV variability, 
which is well modelled as warm Comptonisation in NGC 5548. 
A purely thermal disk blackbody spectrum however is not viable 
as it does not reach the soft X-rays. Furthermore, the correlated 
variability of the soft excess component with the UV disk emis¬ 
sion suggests the soft excess is ‘real’ emission and not an arte¬ 
fact of any relativistically-blurred absorption. Such X-ray ab¬ 
sorbing winds are largely transparent at wavelengths longer than 
the Fyman limit, and thus cannot contribute to the optical/UV/X- 
ray correlated variability. In this work, all the X-ray absorption 
has been fully taken into account based on RGS spectroscopy 
of their features. In a recent study by Done et al. (2012), the 
warm Comptonisation origin for the soft excess has been tested 
by further modelling. The authors built an energetically self con¬ 
sistent accretion flow model: an optic ally-thick accretion disk 
is present down to a radius rcorona, below which the accretion 
power is shared between an optically-thick warm corona and an 
optically-thin hot one, with both coronae Compton up-scattering 
the disk emission. In this proposed picture, our fit of the entire 
NIR/optical/UV continuum and the soft excess, without the need 
for an additional thermal component from the disk, suggests that 
the warm corona completely cover the optically-thick accretion 
disk, with a geometry similar to the one proposed for Mrk 509 by 
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Petrucci et al. (2013). We emphasise the importance and useful¬ 
ness of simultaneous multi-wavelength data (like those obtained 
in Mrk 509 and NGC 5548 campaigns) in investigating the na¬ 
ture of the soft excess, as well as other intrinsic components of 
AGN spectra. 

8.2. The continuum at unobscured and obscured epochs 

The basic nature of the underlying broadband continuum in both 
unobscured and obscured epochs is unchanged. However, the 
value of r of the power-law in the stacked 2013 spectrum is 
about 1.61, which is lower than those from archival 2000 and 
2001 ranging over 1.81-1.88. We note this spectral hardening 
observed in 2013 is not an artefact caused by the spectral fit¬ 
ting of the complex absorption by the obscurer. The hard X-ray 
data obtained by NuSTAR and INTEGRAL confirm this spectral 
hardening (Ursini et al. submitted). Additionally, Swift data from 
unobscured epochs of 2005 and 2007 show that T has previously 
been between 1.5-1.6 (Mehdipour et al. in prep). The long-term 
variability of T, presented in Mehdipour et al. (in prep), shows 
that even at some obscured epochs (2012 and 2014) T was con¬ 
sistently higher, ranging typically between 1.7-1.9, whereas dur¬ 
ing the time of our summer 2013 XMM-Newton campaign it was 
lower at 1.5-1.6. So the low T of 1.61 found for the stacked 2013 
spectrum is not related to the fact that the source is obscured and 
that the full effect of its absorption was not correctly taken into 
account. 

8.3. Effect of obscuration on ionisation and thermal stability 
of photoionised gas 

Determining the shape of the ionising SED using a proper 
broadband continuum model is important in assessing the ion¬ 
isation balance of photoionised gas. Eor instance, Chakravorty 
et al. (2012) have shown that just the temperature of the 
accretion disk and strength of the soft X-ray excess component 
can significantly influence the structure and stability of the 
warm absorbers in AGN. In Sect. 7 we presented the cooling 
curves of photoionised gas in NGC 5548 based on broadband 
ionising SEDs of the source for unobscured and obscured cases. 
On the cooling curves, the heating rate is equal to the cooling 
rate so the gas is in thermal equilibrium. This is reached when 
photoionisation and Compton heating are balanced by recom¬ 
bination and Compton cooling, amongst other less dominant 
heating/cooling mechanisms, which are all taken into account 
by the CLOUDY code. As shown in Eig 13, the cooling curves 
from unobscured 2000 and 2001 SEDs display similar shapes. 
But the 2013 curves and in particular the obscured version 
deviate strongly from those corresponding to the unobscured 
SEDs. 

All five cooling curves in Eig 13 show that up to ^ of about 
200 (i.e. \ogf ~ 2.3), ions would fall on a stable branch of the 
cooling curves. Eor Ee specifically, this means up to and includ¬ 
ing Ee xvii ions are stable on the cooling curves corresponding 
to the four unobscured SEDs. But for the cooling curve of the 
obscured SED, only up to Ee viii ions fall on the stable branch of 
the cooling curve. Compared to cooling curves of unobscured 
SEDs, the obscured version in 2013 has a remarkably extended 
unstable branch, where the slope is negative, and so the unstable 
branch covers a broader range in S. Moreover, below a tempera¬ 
ture of 0.1 keV, the obscured cooling curves are shifted to higher 
S and above 0.1 keV to lower S, illustrating a more pronounced 


S-curve. This contrasting shape of the cooling curve is mostly 
induced by the relative lack of ionising EUV/soft X-ray photons 
in the obscured case, changing the interplay between the heating 
and cooling processes. Eurthermore, there are some deviations 
between the 2013 unobscured cooling curve and those from 
2000 and 2001. The main reason for this is spectral hardening 
of the primary power-law component in 2013 (T = 1.6) relative 
to 2000 and 2001 (T =1.8-1.9). This creates more emission at 
hard X-rays affecting the interplay between Compton heating 
and cooling. 

Interestingly, in all five versions of the stability curves dis¬ 
played in Eig. 13, the peak concentration of highly-ionised ions 
of Ee XXIII to Ee xxvi fall on an unstable branch. This means the 
existence of these ions is expected to be significantly diminished 
in photoionised gas of NGC 5548, regardless of obscuration ef¬ 
fects. It is worth noting that the ionising SEDs used for the com¬ 
putation of the cooling curves are based on the direct accretion- 
powered emission from the disk as obtained from the broadband 
continuum modelling. So although the continuum model of the 
SED includes IR emission from the disk, it does not include re¬ 
processed IR emission from the AGN dusty torus which is lo¬ 
cated further out. Therefore, enhanced cooling through inverse 
Comptonisation of the torus IR photons is not included. The ob¬ 
scurer in NGC 5548 is located at distances of only a few light 
days from the central source (Kaastra et al. 2014), so it is ex¬ 
pected to receive only direct disk radiation. On the other hand, 
the more distant warm absorbers at pc scales in NGC 5548 (Arav 
et al. 2014) might be exposed to IR radiation from the torus, 
thus influencing their cooling, which is not taken into account 
in the calculation of the cooling curves here. The cooling curves 
can be used to investigate whether multiple ionisation phases of 
the warm absorber can co-exist in pressure equilibrium or not, 
since phases with overlapping values of E on the curves would 
be in pressure equilibrium. In a forthcoming paper on our cam¬ 
paign (Ebrero et al. in prep) the analysis of different phases of 
the warm absorber and their long-term variability and position 
on the relevant cooling curves are investigated. 


9. Conclusions 

In this work we have determined the intrinsic continuum of 
NGC 5548 from near-infrared (8060 A) to hard X-ray (200 keV) 
energies, using stacked simultaneous XMM-Newton, NuSTAR, 
Swift, INTEGRAL, HST COS, WO and OCA observations from 
our summer 2013 campaign. We have also determined the un¬ 
derlying broadband continuum from XMM-Newton observations 
taken in 2000 and 2001 before NGC 5548 became obscured, 
which helped us in establishing the proper model for the contin¬ 
uum as the X-ray spectrum during the obscured epoch is more 
suppressed and complex. We have taken into account various 
NIR/optical/UV/X-ray emission and absorption processes taking 
place in our line of sight towards the central source, including the 
effects of X-ray obscuration and warm absorption, to uncover the 
intrinsic continuum components. We conclude that; 

1. The soft X-ray excess in NGC 5548 can be explained 
as the tail of a continuum component extending from 
NIR/optical/UV to soft X-ray energies, produced by 
Compton up-scattering of the thermal seed photons from the 
accretion disk, in a warm (Tg 0.15 keV), optically thick 
(t 23) corona as part of the inner disk. 
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2. The NIR/optical/UV continuum is consistent with being 
composed of a single Comptonised disk component which 
also produces the soft X-ray excess. There is no evidence 
of an additional purely-thermal disk component or an 
additional reprocessed component from the disk. 

3. The nature of the underlying broadband continuum in both 
unobscured and obscured epochs is unchanged. However, 
in 2013 the primary power-law displays a harder spectrum 
(T = 1.6) than in 2000 and 2001, but this has no apparent 
relation to the obscuration event. 

4. The heavy absorption of the ionising SED EUV/soft X-ray 
photons by the obscurer has a significant impact on the 
ionisation balance of the more distant photoionised gas 
(warm absorbers), producing an extended thermally unsta¬ 
ble branch on the cooling curve. 

5. The thermal stability analysis of photoionised gas in 
NGC 5548 indicates highly-ionised ions of Eexxiii to 
Ee XXVI fall on an unstable branch of cooling curves, regard¬ 
less of the obscuration effects. 
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Appendix A: Observations and data reduction 

A.1. XMM-Newton 
EPIC 

The European Photon Imaging Camera (EPIC - Striider et al. 
2001 ; Turner et al. 2001 ) data of NGC 5548 were first processed 
using the Science Analysis System (SAS vl3.0). The EPIC 
instruments were operating in the Small-Window mode with 
the thin-filter applied. Periods of high flaring background were 
filtered out by applying the #XMMEA_EP and #XMMEA_EM filter¬ 
ings for pn and MOS, respectively. This is done by extracting 
single event, high energy light curves, in order to create a 
set of good-time-intervals (GTIs) for use in conjunction with 
the internal GTI tables, to exclude intervals of flaring particle 
background exceeding 0.4 count/s for pn and 0.35 count/s for 
MOS. The EPIC spectra and lightcurves were extracted from 
a circular region centred on the source with a radius of 40". 
Eor pn the background was extracted from a nearby source-free 
region of radius 50" on the same CCD as the source, whereas 
for MOS a background radius of 100" was used from a CCD 
other than that of the source. Eor archival (2000, 2001) EPIC 
data, there was weak pile-up (few % for pn); however, for the 
2013 observations due to the suppressed X-ray fluxes no pile-up 
was detected. The single and double events were selected for 
both the pn (PATTERN <= 4) and MOS (PATTERN <= 12). The 
EPIC lightcurves were background-subtracted and corrected 
(using the epiclccorr SAS task) for various effects on the 
detection efficiency. Response matrices were generated for the 
spectrum of each observation using the rmfgen and arfgen 
tasks. With EPIC-pn operating in the Small-Window mode, 
the final cleaned EPIC-pn exposure for each XMM-Newton 
observation was on average about 37 ks, which is about 67% of 
the duration of an observation (i.e. about 55 ks as given in Table 
1 ). 

The 2013-2014 EPIC-pn data of NGC 5548 were found to be 
affected by larger than expected long-term degradation/evolution 
of the charge transfer inefficiency (CTI). This creates a shift of 
about 50-60 eV at 6 keV. The CTI problem had previously been 
found in the Mrk 509 campaign EPIC-pn data (Ponti et al. 2013), 
and so the long-term CTI was corrected and implemented in new 
SAS CCEs. Eor the new EPIC-pn data of NGC 5548, further 
CTI corrections had to be applied and the use of the new CCEs 


were not sufficient to correct for all remaining CTI. We thus used 
an ad-hoc correcting energy shift as described in more detail in 
Cappi et al. (in prep), in which the full data analysis and mod¬ 
elling of the EPIC data is reported. This energy shift results in 
a poor fit near the energy of the gold M-edge of the telescope 
mirror, therefore, the 2.0-2.4 keV region was omitted from our 
spectral modelling of the 2013 EPIC-pn data. 

RGS 

The RGS (den Herder et al. 2001) instruments were operated 
in the standard Spectro-)-Q mode for the XMM-Newton obser¬ 
vations of NGC 5548. The processing of the RGS data has 
been performed at a more advanced level than the standard SAS 
pipelines. The details of this enhanced RGS data reduction tech¬ 
nique are reported in Kaastra et al. (2011b). The procedure 
takes advantage of the multi-pointing mode of XMM-Newton 
and utilises accurate relative calibration for the effective area of 
the RGS. It also incorporates an accurate absolute wavelength 
calibration to improve the method for stacking time-dependent 
RGS spectra and enhancing the efficiency of the spectral fitting. 
The Kaastra et al. (201 lb) procedure was used to produce a com¬ 
bined RGS-l/RGS-2 spectrum and response matrix (in SPEX for¬ 
mat) for each XMM-Newton observation of NGC 5548, includ¬ 
ing archival 2000 and 2001 observations. Also using this proce¬ 
dure a single stacked RGS spectrum and response matrix from 
the 12 XMM-Newton observations of summer 2013 was created 
for spectral fitting. 

OM 

The Optical Monitor (OM - Mason et al. 2001) photometric 
filters were operated in the Science User Defined image/fast 
mode for the 2013-2014 observations and in the Image mode 
for the archival observations of NGC 5548. In each 2013-2014 
observation, OM images were taken with the V, B, U, UVWl, 
UVM2, and UVW2 filters, with an exposure time of 4.4 ks for 
each image. The OM images of NGC 5548 were processed with 
SAS vl3.0 omi chain pipeline with the standard parameters, 
to apply the necessary corrections including the removal of 
Modulo-8 fixed pattern noise, and produce images ready for 
aperture photometry. The aperture photometry on each image 
was done using the omsource program. Eor those archival 
datasets in which exposures in some of the six OM filters have 
not been taken, we used the relations derived in Mehdipour et 
al. (in prep) between different photometric filters from extensive 
Swift monitoring of NGC 5548, to calculate the missing fluxes 
in order to carry out a uniform and consistent spectral fitting 
of the data. Source and background regions were selected to 
extract count rates, with the necessary corrections applied, 
including those for the point spread function, coincidence 
losses and time-dependent sensitivity. The OM count rates were 
extracted from a circle of 5.0" centred on the source nucleus. 
The background was extracted from a source-free region of 
the same size. The OM count rates were then transformed into 
the standard OGIP EITS format (Arnaud et al. 1992) to be 
used together with their corresponding response files in X-ray 
spectral fitting packages. Lastly, the files were converted into 
the SPEX format using the auxiliary trafo program of the SPEX 
package. 

In addition to the OM photometric filters, spectra with the 
Optical and UV OM grisms were also obtained in each 2013 
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XMM-Newton observation in the Full Low mode. The exposure 
time of each grism spectrum is 5.0 ks. The grism data were 
first processed with the SAS vl3.0 omgchainpipeline. The nec¬ 
essary corrections, including that for Modulo-8 fixed pattern 
noise and removal of scattered light features, were applied to 
obtain undistorted and rotated grism images. We then used the 
omg source program to interactively identify zero and first order 
dispersion spectra of our source, properly define the source and 
background extraction regions and extract the calibrated spec¬ 
trum from the grism images of each observation. For the Optical 
grism the extraction and calibration was extended beyond the 
standard results to just below 7000 A as described in Mehdipour 
et al. (2011) in order to capture the Ha /16563 line. Similar to the 
OM broadband filter data, the grism spectra were also converted 
into SPEX format for simultaneous spectral modelling with other 
data. 

A.2. Swift 
XRT 

For the Swift X-ray Telescope (XRT - Burrows et al. 2005) data 
reduction, we used the procedure detailed in Evans et al. (2009), 
which is an enhanced version of the standard Swift processing in¬ 
cluding several modifications. This tool is made available online 
on the UK Swift Science Data Centre (UKSSDC ^). The XRT in¬ 
strument operated in the Photon Counting (PC) mode for nearly 
all of the NGC 5548 observations. The XRT data have been cor¬ 
rected forbad pixels and effects of vignetting and PSF to produce 
cleaned event files using the xrtpipeline script. The data were 
not significantly affected by pile-up at the observed low count 
rates of the source. Using the event lists, exposure maps were 
made, which together with the source spectrum, were used to 
create a corresponding Ancillary Response File (ARF) for each 
interval. The corrections have been applied to both spectra and 
lightcurves. The optimum extraction radius for data products de¬ 
pends on the PC-mode count rate as reported in Evans et al. 
(2009), which for NGC 5548 translated most of the time into 
an extraction radius of 25 pixels (59.0"). The XRT lightcurves 
at different energies were constructed from each Swift snapshot 
observation of NGC 5548. The default grades of 0-12 in the 
PC mode were used for event selections. The corresponding re¬ 
sponse matrix files (RME) for the Swift observations were ob¬ 
tained from the HEASARC Calibration Database (CALDB). 

UVOT 

The Swift UV/Optical Telescope (UVOT - Roming et al. 2005) 
data of NGC 5548 from Image-mode operations were taken 
with the six primary photometric filters of V, B, U, UVWl, 
UVM2 and UVW2. The uvotsource tool was used to perform 
aperture photometry using a circular aperture radius of 5.0". The 
standard instrumental corrections and calibrations according to 
Poole et al. (2008) were applied. Any data from when Swift 
had bad tracking were removed. The continuous monitoring 
of NGC 5548 indicated occasional exposures in UVOT filters 
in which the flux suddenly dropped by about 20% or more 
compared to the adjacent mean; these exposures were found to 
be clustered at particular regions of the detector and their cause 
and fix is currently being investigated by the UVOT calibration 
team. These occasional observations have been Altered out of 
our analysis (about 2% of all UVOT observations). Eor the 
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purpose of spectral fitting with SPEX, the Alter count rates and 
the corresponding response matrices in each Alter were created. 

In addition to the photometric exposures, there was an exten¬ 
sive monitoring of NGC 5548 with the UV grism of UVOT be¬ 
tween 1 April and 12 September 2013. The grism observations 
of NGC 5548, which were taken about every two days, were 
performed in the Clocked mode, which is the optimum mode for 
low background as well as no presence of zeroth orders in the up¬ 
per left half of the detector. The Swift UVOT Grism Calibration 
and Software UVOTPY^ vl (Kuin et al. in prep) was used to pro¬ 
duce the UV grism spectra of NGC 5548. This software is more 
enhanced than the standard tools available within HEASOET 
The uvotgraspcorr program was used to aspect correct the 
UVOT grism images. The uvotgetspec module of UVOTPY was 
used to extract flrst-order spectra of NGC 5548 and also produce 
the corresponding RME response matrices. The grism spectra 
were then converted into SPEX format for spectral modelling. 

A.3. NuSTAR 

The NuSTAR observations were reduced using the NuSTAR Data 
Analysis Software (NUSTARDAS) vl.3.1, which was utilised 
as part of HEASOET v6.14 distribution. Eor instrumental cal¬ 
ibration of our NuSTAR data, CALDB version v20131007 was 
used. The data were processed with the standard pipeline script 
nupipeline to produce level 1 calibrated and level 2 cleaned 
event Ales. The data from the South Atlantic Anomaly passages 
have been Altered out and event Ales were cleaned with the stan¬ 
dard depth correction, which signiflcantly reduces the internal 
background at high energies. The source was extracted from a 
circular region (radius ~ 110"), with the background extracted 
from a source-free area of equal size on the same detector. Then 
the nuproducts script was run to create level 3 products (spec¬ 
tra, lightcurves, ARE and RME response Ales) for each of the 
two hard X-ray telescope modules (EPMA and EPMB) onboard 
NuSTAR. The spectra and corresponding response Ales of the 
two telescopes were combined for spectral modelling using the 
mathpha, addrmf, and addarf tools of the HEASOET package. 

A.4. INTEGRAL 

The INTEGRAL IBIS and ISGRI data were reduced with the 
Off-line Scientiflc Analysis (OSA) vlO.O software. The spec¬ 
tra were extracted using the standard spectral extraction script 
ibis_science_analysis. The total INTEGRAL ISGRI spec¬ 
trum for each observation, and the associated RME response 
and ARE ancillary Ales, were then generated using the OSA 
spe.pick tool. 

A.5. Chandra LETGS 

The three Chandra observations were taken with LETGS using 
the HRC-S camera. The LETGS data were flrst reduced using 
the Chandra Interactive Analysis of Observations (CIAO v4.5) 
software to create level 1.5 event Ales. The rest of the processing, 
until the creation of Anal level 2 event Ales, was carried out using 
an in-house software Arst described by Kaastra et al. (2002). This 
procedure foAows the same steps as in the standard pipelines 
of CIAO, but results in improved wavelength and effective area 
calibration. The outcome of this procedure is the production of 


^ http://www.mssl.ucl.ac.uk/-npitik/Grism 
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spectral files and their corresponding response matrices in the 
SPEX format (Kaastra et al. 1996). 

A.6. HSTCOS 

The six HST COS observations were taken concurrently with 
XMM-Newton Ohs. 1, 4, 8, 11, 12 and 14. Each COS observa¬ 
tion consists of an HST two-orbit visit. The observations through 
the Primary Science Aperture were taken with gratings G130M 
(1770 s per visit) and G160M (2215 s per visit), covering the far- 
UV spectral range from 1132 to 1801 A. In addition to the rou¬ 
tine processing of data with the calibration pipelines of STScI, 
further enhanced calibrations (as described in Kriss et al. 2011 
for Mrk 509) were applied to produce the best-quality COS spec¬ 
tra possible. They include refined flux calibrations that take into 
account up-to-date adjustments to the time-dependent sensitivity 
of COS, improved flat-held corrections and wavelength calibra¬ 
tion and an optimal method for combining exposures comprising 
a single visit. More details on the COS calibration of NGC 5548 
data are given in Kaastra et al. (2014), and in Kriss et al. (in 
prep), in which the full analysis of the HST COS spectra are also 
reported. For the purpose of modelling the UV continuum as part 
of our broadband continuum modelling, we extracted the COS 
flux from six narrow spectral bandpasses, which are free of emis¬ 
sion and absorption lines: 1155-1165, 1364.5-1369.5, 1470- 
1490, 1730-1740, 1730-1770 and 1785-1800 A. In order to 
use these measurements in our simultaneous spectral fitting with 
other optical/UV and X-ray data, we converted the COS fluxes 
into count rate units used in the OGIP spectral file format. To do 
this we used the time-dependent sensitivity curves as a function 
of wavelength, obtained for the COS data in Kriss et al. (in prep), 
which are in units of (count s“' pixeE')/(erg cm“^ s“^ A“*), to 
calculate the count rate and the effective area over the six narrow 
bandpasses. For the G130M and G160M modes, the pixel size is 
about 0.010 A and 0.012 A, respectively. 

A. 7. OCA 

The Observatorio Cerro Armazones (OCA) monitoring of 
NGC 5548 was carried out between 20 May and 25 July 
2013. For our observations the 25 cm Berlin Exoplanet Search 
Telescope II (BEST-II) was utilised, which is equipped with a 
Peltier-cooled 4096 x 4096 pixel Finger Fakes Imager CCD 
KAF-16801. It has a fleld-of-view of 1.7 x 1.7 degrees with a 
pixel size of 9 fim. Photometric measurements were taken in the 
optical Johnson B, V and R Alters. For each Alter, 6 dithered im¬ 
ages were acquired, with 150 s exposure time for each individual 
image. There were observations on 27 days, from which 24 had 
exposures in B, 19 in R, and 7 in the V Alter. Data reduction and 
processing of the OCA data was performed in the same manner 
as reported in Pozo Nunez et al. (2013). 

A.8. Wise Observatory 

The Wise Observatory (WO) monitoring of NGC 5548 was done 
with the Centurion 46-cm telescope, and images were obtained 
with the standard Bessel B, V, R, and I Alters. They cover the 
time from beginning of June to end of September 2013. There 
were also further monitoring from December 2013 to July 2014. 
The observations were done with STF-6303 CCD which gives 
a held of view of 70x50 arcmin and each exposure time was 
300 s. Standard data reduction was carried out with the Image 
Reduction and Analysis Facility (IRAF) packages using PSF 


photometry. Each night a few exposures in each Alter were ob¬ 
tained and after confirming that there were no intra-night varia¬ 
tions within the uncertainties, they were averaged to produce one 
measurement for each night. Calibration to absolute magnitudes 
was done in each Alter using several stars in the held which their 
absolute magnitudes were taken from the SIMBAD database. 
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